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(1)  Foreword 

The  main  goal  of  this  program  was  to  design,  fabricate  and  test  a  semiconductor  device  capable  of 
demonstrating  the  fundamental  physics  required  for  the  realization  of  a  spin-coherent,  single  photon 
transmitter/receiver  system.  These  requirements  include  first,  careful  tailoring  of  the  g  factor  for  conduction 
band  electrons  in  the  InGaAs/InP  materials  system,  and  second,  the  confinement,  detection  and  manipulation  of 
single  electrons  under  lithographically  defined  gates.  In  this  final  report,  we  will  review  our  original  program 
objectives  and  describe  our  major  accomplishments.  In  Section  (4),  we  will  summarize  our  accomplishments  in 
a  timeline  chart,  which  will  clearly  show  how  our  actual  rate  of  progress  compared  with  our  proposed  plan. 
Section  (5)  will  detail  the  major  accomplishments  called  out  in  the  timeline  chart.  Later,  we  will  present  a  self- 
evaluation,  identifying  the  lessons  learned  and  what  we  would  have  done  differently.  Finally,  we  will  conclude 
by  discussing  the  outstanding  challenges  left  for  semiconductor-based  quantum  information  science  and 
technology  and  offer  our  assessment  of  the  future  of  this  field. 
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(4)  Statement  of  the  problem  studied 


In  our  original  proposal,  we  described  a  program  to 
validate  the  basic  physics  of  a  solid-state  quantum 
repeater  based  on  a  design  conceived  at  UCLA  by 
Prof.  Eli  Yablonovitch  (Fig.  1).  The  idea  was  that  it 
should  be  possible  to  build  a  solid-state  device 
consisting  of  simply  a  single-photon  detector,  three 
single-qubit  logic  gates,  and  a  single-photon  emitter, 
capable  of  executing  the  teleportation  algorithm  on  an 
incident  photon,  and  faithfully  transmitting  its 
quantum  state  to  another  repeater  connected  using 
conventional  optical  fiber.  We  argued  that  the  device 
would  work  by  first  converting  the  incident  photon 
polarization  state  into  an  electron  spin  state,  allowing 
it  to  interact  with  the  adjacent  electron  spins  trapped 
under  the  single-qubit  logic  gates,  and  finally  emitting 
an  entangled  photon  from  the  single-photon  emitter. 
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Figure  1:  Schematic  diagram  showing  the  UCLA 
concept  for  the  quantum  repeater  described  in  our 
original  proposal.  The  design  uses  three  single-qubit 
logic  gates  (labeled  ‘A’),  a  single-photon  detector,  and 
a  single-photon  emitter  based  in  the  InGaAs/lnP 
materials  system. 


For  the  program,  we  proposed  to  demonstrate  the  fundamental  physics  of  such  a  device,  (1)  that  the  g  factor  for 
electrons  in  an  InGaAs/lnP  quantum  well  (QW)  can  be  tailored  through  an  appropriate  choice  of  Ga 
composition  and  layer  thickness  and  (2)  that  individual  electrons  can  be  trapped  under  an  electrostatic  gate 
electrode  and  detected  by  monitoring  the  conductance  of  an  integrated  readout  channel.  The  motivation  for 
choosing  these  particular  demonstrations  is  described  below.  First,  the  single-photon  detector  must  be  able  to 
take  an  unknown  quantum  state,  represented  as  a  photon  polarization  and  create  an  electron-hole  pair,  for  which 
the  unknown  state  would  be  transmitted  to  the  spin  state  of  the  electron  (and  not  the  hole).  Therefore,  the  g 
factor  for  electrons  must  be  close  to  zero  and  the  g  factor  for  holes  must  be  large  enough  for  the  hole  to  be 
created  with  definite  spin  [1].  In  a  paper  by  Kiselev  et  al.  [2],  a  process  known  as  g  factor  engineering  is 
described  for  the  case  of  InGaAs  QWs  in  InP.  Through  an  appropriate  choice  of  Ga  composition  and  InGaAs 
layer  thickness,  the  g  factor  of  an  electron  trapped  in  the  InGaAs  layer  was  predicted  to  sample  both  the 
negative  g  factor  of  InGaAs  and  the  positive  g  factor  of  InP,  resulting  in  an  effective  g  factor  of  zero.  We 
therefore  set  as  our  first  objective,  to  demonstrate  g  factor  engineering  in  InGaAs/lnP. 


The  second  objective,  to  design,  fabricate,  and  test  a  single-qubit  logic  gate,  arises  from  the  requirement  that  the 
repeater  concept  requires  three  trapped  electrons,  in  order  to  execute  the  teleportation  algorithm.  We  therefore 
set  out  to  develop  the  theoretical  basis  for  designing  a  top-gated,  electrostatic  quantum  dot  in  InGaAs/lnP.  In 
addition,  HRL  has  considerable  expertise  in  the  growth  and  processing  of  InGaAs/lnP  heterstructures  and  so  we 
proposed  a  complementary  experimental  effort  to  fabricate  and  test  prototype  devices  and  look  for  evidence  of 
single-electron  capture.  The  choice  of  materials  system  was  motivated  by  the  desire  for  the  repeater  design  to 
be  compatible  with  photon  energies  typically  used  in  conventional  optical  networks. 


Below,  in  the  timeline  shown  in  Fig.  2  and  in  bulletized  form,  we  summarize  the  principle  scientific  goals  of  the 
program,  as  laid  out  in  our  original  proposal  and  plot  our  accomplishments  as  they  were  achieved  throughout 
the  course  of  the  program. 


Goals: 

Theoretical: 

•  Program  Year  1  -  2D  Model  for  electron  confinement  in  InGaAs/lnP 

•  Program  Year  2  -  Full  3D  model  for  electron  confinement  and  overlap 
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•  Program  Year  3 

•  Program  Year  4 

•  Program  Year  5 
Experimental 

•  Program  Year  1 

•  Program  Year  2 

•  Program  Year  3 

•  Program  Year  4 

•  Program  Year  5 


Full  3D  model  for  two-electron  confinement 

Full  3D  model  with  band  structure 

Develop  device  simulations  incorporating  uncertainty 

Measure  the  electron  g  factor  in  InGaAs/InP  2D  electron  gases 
Demonstrate  single-electron  sensitivity  in  quantum  dot  devices 
Demonstrate  single-shot  readout 
Demonstrate  spin-dependent,  single-photon  detection 
Demonstrate  spin  exchange  between  two  coupled  dots 


Accomplishments: 

Theoretical 

•  InGaAs/InP  2D  model  of  confinement 

•  Restricted  3D  model  of  top-gated  quantum  dot  device 

•  Full  quantum  mechanical  simulation  of  realistic  top-gated  device  structures  including  readout 
channels 


EXPERIMENTAL 


GOALS 

g  Factor 
Measurement 


Single-Electron 
Device  Demo 


Spin  Measurement 
and  Single-Spin  ESR 


Spin-Dependent 
Single  Photon 
Detection 


Demonstration  of 
Spin-Exchange 


Final  Report 


Mar.  31  2002 

Mar.  31  2003 

Mar.  31  2004  \ 

Mar.  31  2005 

Mar.  31  2006  \ 

Dec.  31  2006  Y 

Program  Year 

Program  Year  2 

Program  Year  3 

Program  Year  4 

Program  Year  5 

Extended  Per. 

2001 


2002 


2003 


2004 


2005 


2006 


Jul 

Jan  Jul 

Jan  Jul 

Jan  Jul 

Jan  Jul  Jan 

iJul 

A 

i  n  o 

_ Q _ 

ono 

i 

i 
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Project 

Start 

July  1, 2001 


InGaAs/InP 
2D  Model  of 
Confinement 

Mar.  31, 2002 


Full  3D  Model  of 
Top-Gated  Device 

Mar.  31, 2003 


Full  3D  Model  for 
2  Electrons 

Mar.  31 , 2004 


Mar.  31, 2005 


Simulations 
Full  3D  Model  \  Including 
,  with  Band  Structure,  Uncertainty 


Mar.  31, 2006 


THEORETICAL 

ACCOMPLISHMENTS 


A 

□ 

O 

o 

o 

o 

o 


InGaAs/InP  2D  Model  of  Confinement 

Restricted  3D  Model  of  Top-Gated  Quantum  Dot  Device 

Fabrication  of  Top-Gated  Air  Bridge  Device  Structures  with  Sub-50  nm  Feature  Sizes 

First  Demonstration  of  g  Factor  Engineered  InGaAs/InP  Quantum  Weils  with  Large  g  Factor  Shifts 

Full  Quantum  Mechanical  Simulation  of  Realistic  Top-Gated  Device  Structures  Including  Readout  Channels 

Demonstration  of  Single-Electron  Sensitivity  in  Top-Gated  InGaAs/InP  Quantum  Dot  Devices 

First  Demonstration  of  Few-Electron  Quantum  Dot  Spectroscopy  using  Random  Telegraph  Analysis 


Figure  2:  Timeline  summarizing  our  original  program  goals  (blue  print)  and  our  accomplishments  (green  markers) 
and  their  approximate  date  of  accomplishment  as  they  were  achieved  throughout  the  course  of  the  program. 
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Experimental 

•  Fabrication  of  top-gated  air-bridged  device  structures  with  sub-50  nm  feature  sizes 

•  First  demonstration  of  g  factor  engineered  InGaAs/InP  QWs  with  large  g  factor  shifts 

•  Demonstration  of  single-electron  sensitivity  in  top-gated  InGaAs/InP  quantum  dot  devices 

•  First  demonstration  of  few-electron  quantum  dot  spectroscopy  using  random  telegraph  analysis 

In  Section  (5),  we  describe  in  more  detail  the  nature  each  of  our  major  accomplishments. 


(5)  Summary  of  the  most  important  results 

(5.1)  Theoretical  accomplishments 

The  prototype  single-qubit  logic  gate  device  we  designed  had  to  satisfy  a  large  number  of  constraints  that  arise 
from  engineering  considerations  as  well  as  basic  physics.  The  device  must  allow  for  the  confinement  of  either 
one  or  zero  electrons  in  an  upper  QW  under  reasonable  gate  voltage  bias  conditions  while  simultaneously 
allowing  for  a  significantly  depleted  (and  nearly  one  dimensional)  2D  electron  gas  (2DEG)  in  the  lower  QW. 
The  first  condition  is  required  in  order  to  form  an  effective  quantum  dot  for  the  electron  whose  spin  state  will 
serve  as  the  qubit.  The  second  condition  arises  from  the  fact  that  the  current  flowing  through  the  lower  QW 
(referred  to  as  the  channel)  is  sensitive  to  the  presence  of  single  charges  in  the  vicinity  of  the  channel  only  when 
severely  depleted  laterally  by  applying  negative  bias  to  a  quantum  point  contact  (QPC).  Since  this  is  our  method 
for  detecting  the  presence  (or  absence)  of  an  electron  in  the  quantum  dot,  satisfying  this  constraint  was  also 
critical. 

(5.1.1)  InGaAs/InP  2D  model  of  confinement 

The  desired  device  has  a  “double  pinchoff,”  by  which  we  mean  that  there  is  a  single  trapped  electron  in  the 
quantum  dot  (in  the  upper  QW)  and  a  single  transverse  state  in  the  quantum  wire  (in  the  lower  QW).  The 
electronic  properties  of  this  system  are  described  by  a  Poisson  equation  for  the  electrostatic  potential  and  a 
Schrodinger  equation  for  the  electron  wave  functions.  Solution  of  the  Schrodinger-Poisson  equations  for  this 
system  provides  a  method  for  assessing  the  double  pinchoff  property.  The  Schrodinger  equation  should  have  a 
single  eigenvalue  below  the  Fermi  energy  in  the  quantum  dot  and  a  single  eigenvalue,  corresponding  to 
transverse  eigenvector,  below  the  Fermi  energy  in  the  QW.  We  used  solutions  of  the  Schrodinger-Poisson 
system  as  a  design  tool  to  find  a  design  that  attains  the  desired  double  pinchoff.  This  effort  was  critical  to 
success  of  the  project,  since  it  would  have  been  considerably  more  difficult  to  find  a  design  by  a  purely 
experimental  approach. 

We  developed  two  approaches  to  solution  of  the  Schrodinger-Poisson  system.  The  first  is  a  direct  numerical 
solution  of  the  equations  using  a  finite  difference  method.  This  has  been  carried  out  in  one  and  two  dimensions. 
The  challenge  of  this  problem  is  that  the  system  is  numerically  stiff  in  the  wells  and  that  the  nonlinear  “self- 
consistent”  terms  (the  contributions  to  the  electrostatic  potential  from  the  electrons  in  the  quantum  dot  and  wire) 
make  the  eigenvalue  problems  nonlinear. 

Our  second  approach  was  an  approximation  based  on  a  combination  of  exact  and  approximate  solutions  of  the 
Poisson  and  Schrodinger  equations,  for  example  the  Schrodinger  equation  with  square  well  and  parabolic 
potentials,  in  one,  two  or  three  dimensions.  This  has  not  included  self-consistent  terms.  From  these 
approximations,  a  somewhat  complicated  but  completely  analytic  fonnula  (including  a  few  iterative  root 
finders)  for  the  eigenvalues  was  obtained. 
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This  latter  method  has  been  implemented  in  a  Matlab  program  that  allows  for  fast  computation  of  proposed 
designs.  As  a  first  step,  we  searched  a  large  design  space  and  found  10  successful  designs  out  of  about  10 
million  trials.  The  resulting  designs  showed  some  surprising  features;  for  example,  the  width  of  the  two  QWs 
were  about  equal,  whereas  we  had  expected  the  lower  well  (the  channel)  to  be  much  wider.  Based  on  these 
results,  we  narrowed  the  search  and  implemented  a  criterion  for  optimality  of  the  design. 

The  direct  numerical  solution  has  been  used  to  validate  the  analytic  model  and  estimate  the  uncertainty  in  the 
results  associated  with  the  analytical  modeling  assumptions.  In  addition,  the  direct  numerical  simulation  was 
used  to  understand  the  role  of  the  choice  of  boundary  conditions  for  the  Schrodinger-Poisson  equation.  These 
boundary  conditions  are  not  clear  from  the  physics,  and  we  typically  used  conditions  of  prescribed  potential  on 
the  top  surface  (corresponding  to  Schottky  barrier  offsets)  and  no  flux  on  the  bottom  surface  (corresponding  to 
zero  electric  field  at  infinity).  Alternate,  seemingly  reasonable  boundary  conditions  exist,  and  the  utility  of  these 
boundary  conditions  needs  to  be  detennined. 


(5.1.2)  Restricted  3D  model  of  top-gated  quantum  dot  device 


The  UCLA  efforts  on  the  direct  numerical  simulation  of  the  coupled  Schrodinger  -Poisson  equations  have 
concentrated  on  the  extension  of  the  one  and  two-dimensional  simulations  to  full  three-dimensional  simulations. 
The  three  dimensional  simulation  of  the  double  well  quantum  device  was  completed  and  was  validated  against 
the  Nextnano3  simulation  tool  that  we  obtained 
in  late  2002  from  the  Theoretical  Semiconductor 
Group  at  the  Schottky  Institute  in  Munich  (Peter 
Vogl,  director).  The  high  computational 
efficiency  of  the  UCLA  simulation  was  obtained 
by  exploiting  the  approximate  separability  of  the 
model  equations.  Also,  during  the  construction  of 
this  simulation,  a  new,  robust,  gradient  free, 
iterative  method  for  solving  the  consistent 
Schrodinger-Poisson  equations  was  developed. 


XtOt" 


dzl 


The  starting  point  for  the  3D  device  design  was  a 
structure  identified  using  the  Caflisch  semi- 
analytic  model.  The  challenge  for  this  device  is 
to  achieve  simultaneous  depletion  of  the  lower 
QW  (as  defined  by  the  presence  of  only  a  few 
longitudinal  quantum-wire-like  states)  and 
occupation  of  the  upper  well  with  a  single  state 
under  low  (<1V)  forward  (positive)  bias 
conditions.  The  semi-analytic  calculations 
produced  optimal  designs  that  were  all  very 
similar;  variations  from  one  optimal  design  to  the 
next  were,  for  the  most  part,  almost  within  the 
experimental  uncertainty  in  the  parameters.  A 
typical  structure  was  input  into  NextNano3  and 
studied  extensively  over  a  range  of  gate  voltages. 
Many  numerical  and  modeling  issues,  such  as 
code  convergence,  were  addressed  during  the 
course  of  this  analysis  and  extensive  feedback 
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Figure  3:  Device  geometry  for  NextNano3  simulation  of  our 
device  structure. 
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was  given  to  the  Vogl  group  related  to  bugs  and  desired  improvements.  The  final  results  from  NextNano3 
indicated  that,  generally,  states  in  the  dot  were  harder  to  achieve,  i.e.  more  forward  bias  on  the  dot  gate  was 
required  for  a  given  QPC  voltage.  This  difference  was  attributed  to  the  approximations  inherent  in  the  semi- 
analytic  model. 

Fig.  3  shows  the  device  geometry  that  was  used  in  the  Nextnano3  simulation.  Periodic  boundary  conditions  are 
assumed  in  the  x-direction,  which  is  equivalent  to  making  the  assumption  that  the  gates  are  infinitely  long.  For 
most  of  the  devices  fabricated,  this  is  not  a  good  assumption  and  the  effect  of  realistic  gate  geometries  was 
studied  later.  Realistic  gate  geometries  can  be  generated  with  Nextnano3,  but  the  large  simulation  size  that  is 
required  makes  the  run  times  impractical  (days  to  converge  a  single  result).  These  simulations  are  possible  with 
the  UCLA  code,  however,  because  it  is  many  times  faster  than  Nextnano3.  In  the  meantime,  we  felt  that  results 
using  this  geometry  would  still  be  perfectly  useful  as  a  first  pass  design  tool.  One  issue  regarding  convergence 
of  the  results,  we  discovered,  is  the  depth  into  the  substrate  below  the  last  delta-doped  layer  (dz7  in  Fig.  3). 
Using  a  background  doping  level  for  InP  of  10'  cm'  (n-type)  suggested  by  our  device  growers,  the  depth  into 
the  substrate  needs  to  be  quite  substantial  before  the  conduction  band  reaches  its  equilibrium  value.  Because 
Nextnano3  allows  variable  grid  spacings  within  a  region,  however,  this  did  not  add  any  significant 
computational  time  to  the  simulations  since  the  conduction  band  is  changing  smoothly  in  that  region  and  a  large 
grid  spacing  can  be  used. 


Figure  4:  Conduction  band  edge  z-  (left)  and  y-dependence  (right)  at  the  x  midpoint  for  the  bias  conditions  given  in  the 
text. 

Fig.  4  shows  the  conduction  band  edge  produced  by  the  simulation  for  a  QPC  voltage  of  -0.53V  and  a  gate 
voltage  of  +0.25V.  The  left-most  figure  is  the  z-dependence  of  the  conduction  band  at  the  x  midpoint  at  various 
y  values  (see  Fig.  3  for  orientation  of  the  axes).  The  right-most  figure  is  the  y-dependence  of  the  conduction 
band,  also  at  the  x  midpoint,  at  different  values  of  z  or  depth  into  the  structure.  The  most  relevant  feature  in  the 
z-dependence  data  is  that  both  wells  are  only  slightly  below  the  Fermi  level  for  y  values  that  correspond  to 
being  nearly  in  between  the  QPC  gates.  The  y-dependence  data  for  different  z  values  is  more  interesting.  Here, 
the  lower  two  curves  are  the  conduction  band  in  the  center  of  the  two  QWs  (upper  well  -  green;  lower  well  - 
purple).  It  is  clearer  here  that  the  conduction  band  in  both  wells  is  below  the  Fermi  level  only  in  between  the 
QPC  gates  and  is  at  a  depth  where  there  are  very  few  states.  Also  noticeable  is  that  the  effect  of  the  dot  gate  is 
very  small  in  both  wells.  On  this  voltage  scale,  it  appears  to  be  a  concern  that  there  is  no  real  confinement  due 
to  the  dot.  On  a  finer  scale,  however,  there  is  a  clear  potential  well  with  a  confinement  depth  of  a  few  meV, 
more  than  deep  enough  to  capture  electrons  at  He3  temperatures  (250  ink). 
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Fig.  5  shows  wavefunction  amplitudes  calculated  by  Nextnano  for  the  bias  conditions  described  above  for 
several  of  the  states  just  below  the  Fermi  level.  The  amplitudes  of  the  wavefunctions  are  color  coded  with  blue 
representing  zero  and  red  representing  the  maximum  value.  They  are  shown  in  the  y-z  plane  at  the  x  midpoint 
with  vertical  being  the  y-direction  and  horizontal,  the  z-direction.  The  left  side  of  each  plot  is  the  top  of  the 
device  structure  and  moving  to  the  right  corresponds  to  moving  down  into  the  structure.  Wavefunctions  that 
have  significant  amplitudes  toward  the  right  side  of  each  plot  are  localized  in  the  lower  QW,  or  channel.  The 
one  plot  (105)  that  has  most  of  the  amplitude  on  the  left  side  represents  the  one  state  of  all  those  below  the 
Fermi  level  that  is  localized  in  the  upper  well.  This  state  represents  two  electrons  (this  calculation  assumes 
degenerate  spin  states  and  no  electron-electron  interaction)  in  the  quantum  dot.  Note  that  there  are  up  to  12 
occupied  states  in  the  y-direction  the  lower  well  for  this  case;  the  reason  there  are  more  than  100  occupied  states 
total  is  due  to  the  transverse  states  in  the  lower  well  that  are  really  extended  states  in  the  x-direction 


91  92  93  94  95  96  97  98  99  100 


Figure  5:  Wavefunction  amplitudes  in  the  y-z  plane  at  the  x  midpoint  for  the  bias  conditions  described  in  the  text. 


The  semi-analytic  calculation  showed  that  the  first  dot  state  was  present  with  fewer  channel  states.  As 
mentioned  previously,  this  is  almost  certainly  due  to  the  approximations  inherent  in  that  calculation.  We  have 
some  concern  that  there  may  to  too  many  states  in  the  channel  before  the  first  dot  state  is  present  -  the  presence 
of  many  states  in  the  lower  well  indicates  that  the  channel  is  not  close  to  depletion.  It  remains  an  open  question 
whether,  for  this  device,  the  channel  needs  to  be  almost  completely  depleted  (i.e.  only  one  or  two  states  present) 
in  order  to  be  sensitive  to  single  electrons.  The  channel  is  only  lOnm  from  the  electron,  so  we  believe  that  it  will 
be  sensitive  to  single  charges  even  when  the  channel  is  not  near  full  depletion. 
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Figure  6:  Quantum  well  occupation  as  a  function  of  assumed  Schottky  barrier  and  effective  doping  activation 


In  order  to  parameterize  our  model  with  the  experimental  results,  we  calculated  the  QW  occupations  of  several 
one-dimensional  systems  with  Nextnano3,  at  a  number  of  different  doping  values,  assuming  100%  activation  of 
the  dopants.  We  also  allowed  for  a  wide  range  of  surface  potentials,  which  could  be  interpreted  as  Schottky 
barriers  plus  a  uniform  applied  gate  voltage.  The  results  are  shown  in  Fig.  6.  In  this  figure,  the  lower  well 
occupation  is  indicated  by  the  red  color  level,  and  the  upper  well  by  the  green  color  level.  When  both  wells  are 
occupied,  the  figure  shows  a  yellow  color.  The  contour  lines  are  shown  separately  for  the  lower  well  (red)  and 
the  upper  well  (green).  The  presence  of  lower  well  charge  only  indicates  being  up  in  the  bright  red  region  of  the 
graph  in  Fig.  6,  and  is  consistent  with  being  somewhere  between  about  50  and  no  more  than  200  on  the 
horizontal  axis.  Our  presumption  that  the  Schottky  value  was  around  0.75  V  would  imply  that  the  charge  density 
in  the  lower  well  was  very  near  the  nominally  designed  value. 


(5.1.3)  Full  quantum  mechanical  simulation  of  realistic  top-gated  device  structures  including  readout 
channels 

At  the  end  of  Program  Year  4,  we  completed  the  implementation  of  a  fully  three-dimensional  Poisson- 
Schrodinger  solver  for  planar  heterostructures  with  patterned  metallic  surface  contacts.  As  distinct  from  the 
code  developed  at  UCLA,  this  code  is  the  baseline  code  into  which  we  incorporated  electron-electron 
interactions  that  allowed  us  to  study  in  detail  the  spectroscopy  of  our  gated  quantum  dots. 

The  self-consistent  electrostatic  potential  and  charge  densities  were  computed  following  the  method  of  Trellakis 
et.al.  [3],  in  which  a  predictor-corrector  scheme  is  employed.  Specifically,  the  electrostatic  potential  is  solved 
via  the  Poisson  equation  in  which  a  trial  (predictor)  source  charge  density  is  used.  The  actual  charge  density  is 
then  computed  via  the  Schrodinger  equation  (corrector)  given  this  electrostatic  potential.  The  predictor  is 
derived  from  the  fully  quantum  mechanical  charge  density  by  introducing  a  variation  of  the  local  Fermi  energy 
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proportional  to  the  local  value  of  the  current  electrostatic  potential.  This  predictor  is  thus  a  non-linear  functional 
of  the  electrostatic  potential  and  the  resulting  non-linear  Poisson  equation  is  solved  via  Newton’s  method.  The 
electrostatic  potential  obtained  is  then  used  to  compute  a  new  charge  density  via  the  Schrodinger  equation,  the 
‘corrector’  in  the  overall  scheme.  The  eigenvalue  problem  associated  with  the  discretized  Schrodinger  equation 
is  the  single  most  time  consuming  part  of  this  scheme.  We  address  this  problem  by  first  solving  for  the  charge 
density  quantum  mechanically  only  in  regions  of  our  geometry  for  which  discrete  quantum  states  are  of  interest, 
namely  the  QW/dot  regions.  Second,  we  use  an  eigenvalue  solver  (ARPACK)  which  performs  an  iterative 
eigenvalue  solve  allowing  us  to  only  extract  the  lower  lying  eigenvalues  of  interest.  We  also  worked  with  Chris 
Anderson  at  UCLA  to  incorporate  his  more  efficient  eigenvalue  solver  tailored  specifically  for  this  problem. 
We  have  implemented  an  exchange-correlation  potential  in  the  local  density  approximation  (LDA)  as  a  first 
step  towards  incorporating  more  exact  many-body  effects  in  this  code.  This  is  added  as  an  additional  local 
source  tenn  in  the  Schrodinger  equation. 

A  unique  aspect  of  our  code  is  the  discretization  employed.  The  Schrodinger  and  Poisson  equations  are 
discretized  using  a  high-order  volume  finite  element  scheme.  Arbitrary  order  finite  elements  and  solution 
accuracies  are  selectable  via  specially  constructed  numerical  basis  functions  and  quadrature  rules.  This  will 
allow  for  more  accurate  solutions  to  be  obtained  with  fewer  unknowns,  a  serious  concern  when  solving  three- 
dimensional  problems. 

In  previous  reports,  we  highlighted  the  fact  that  the  major  computational  bottleneck  in  the  three  dimensional 
simulation  is  the  computation  of  the  3D  quantum  states  and  their  associated  energies  (e.g.  the  eigenvectors  and 
eigenvalues  of  the  3D  Schrodinger  operator).  Most  of  the  numerical  work  performed  recently  was  dedicated  to 
reducing  this  computational  bottleneck.  Several  advances  were  made,  the  most  notable  were 

•  The  construction  of  a  robust  iterative  eigevalue/eigenvector  method;  a  method  that  is  capable  of  handling 
problems  in  which  the  eigenvalues  are  very  closely  spaced  or  have  multiplicity  greater  than  one. 

•  The  development  and  implementation  of  an  adaptive  basis  technique  to  accelerate  the  computation  of  states 
associated  with  separable  approximations  to  the  Schrodinger  operator. 

•  The  development  and  implementation  of  an  adaptive  basis  technique  to  accelerate  the  computation  of  states 
of  the  non-separable  Schroedinger  operator. 


The  development  of  the  iterative  method  was  necessitated  by  the  fact  that  the  states  that  must  be  computed  in 
the  simulation  can  have  energies  with  multiplicity  greater  than  one  and  are  very  close  to  each  other  in 
magnitude.  The  multiplicity  arises  because  of  the  symmetry  of  the  geometry  and  the  close  spacing  arises 
because  we  are  approximating  a  continuous  spectrum  when  computing  the  states  contributing  to  the  quantum 
wire  in  the  lower  well.  Both  of  these  features  cause  problems  when  standard  iterative  methods  are  used  and  so 
we  developed  a  new  method  by  combining  two  existing  procedures;  subspace  iteration  with  polynomial 
filtering.  The  method  has  the  advantage  of  simplicity  —  it  is  much  less  complicated  than  Lancoz’s  method  and 
it  does  not  require  any  form  of  “pre-conditioning”  as  Davidson’s  method  or  Inverse  Subspace  Iteration  requires. 
Moreover,  the  method  is  implemented  as  a  C++  template  requiring  only  the  specification  of  an  operator  class 
and  a  vector  class.  This  form  of  C++  packaging  has  enabled  us  to  use  quickly  incorporate  the  procedure  into  all 
of  the  existing  simulations.  The  report  describing  the  procedure  is  in  the  process  of  being  written  up. 

The  development  of  a  robust  iterative  eigenvalue/eigenvector  solver  goes  a  long  way  towards  removing  the 
computational  bottleneck,  but  computing  large  numbers  of  states  associated  with  systems  that  have  hundreds  of 
thousands  of  unknowns  is  still  very  computationally  demanding.  One  way  to  reduce  the  operation  count  is  to 


10 


use  a  basis  that  can  approximate  the  desired  states,  but  has  a  dimension  much  smaller  than  the  space  associated 
with  the  standard  basis.  For  example,  if  we  knew  the  M  states  with  energies  less  than  the  Fermi  energy,  and 
used  those  states  as  a  basis,  we  could  reduce  the  problem  to  one  of  dimension  M.  If  M  is  on  the  order  of 
hundreds,  or  thousands,  the  eigenvalue/eigenvector  problem  can  be  solved  very  quickly.  Unfortunately,  we 
don’t  know  the  states  with  the  lowest  energies  —  that  is  what  we  are  trying  to  find  —  but  we  do  know  the 
potential.  As  a  first  approximation,  the  potential  can  be  used  to  detennine  the  support  of  the  wave  functions 
with  lowest  energies,  and  so  the  basis  we  use  for  the  eigensystem  computation  is  one  whose  support  is 
determined  by  the  potential.  In  the  devices  being  investigated  there  is  strong  confinement  in  the  vertical 
direction  and  in  one  transverse  direction,  so  the  basis  set  used  has  a  much  smaller  dimension  than  the  complete 
system  and  the  computation  time  for  the  eigenvector/eigenvalue  system  is  considerably  reduced. 

The  above  advances  in  the  state  computation  have  resulted  in  a  computationally  efficient  simulation  that 
provides  design  quality  information  for  the  quantum  bit  devices  under  investigation.  Sample  results  that 
demonstrate  the  operational  behavior  of  the  device  are  depicted  in  Figs.  7  and  8.  In  Fig. 7,  an  isosurface  of  the 
electron  density  is  embedded  in  a  visualization  of  the  potential  of  the  sample  device.  In  this  figure,  the  side 
gates  are  biased  so  that  there  is  a  quantum  wire  in  the  lower  well  and  no  states  (e.g.  dot)  in  the  upper  well.  In 
Fig.  8,  we  show  the  electron  density  embedded  in  the  device  with  the  center  gate  biased  to  induce  a  dot  in  the 
upper  well.  (The  isosurface  of  the  dot  in  the  upper  well  is  of  the  same  value  as  the  lower  well,  but  it  is  colored 
red  to  highlight  its  appearance.). 


Figure  7:  The  potential  with  the  quantum  wire  Figure  8:  The  potential  with  the  quantum  wire  and 

embedded  as  an  isosurface.  quantum  dot  embedded  as  isosurfaces. 

(5.2)  Experimental  accomplishments 

Throughout  the  course  of  this  program,  we  have  focused  on  two  principal  experimental  objectives.  The  first,  to 
measure  the  g  factor  of  electrons  confined  to  InGaAs  QWs  in  InP,  depended  critically  on  our  ability  to  model 
and  fabricate  devices  (gated  Hall  bars)  exhibiting  a  high  carrier  mobility  and  a  controllable  electron 
concentration  (through  the  application  of  a  gate  voltage  or  by  adjusting  the  dopant  concentration).  Our 
approach  was  to  induce  a  resonance  between  the  spin-split  Landau  levels  of  a  2DEG  in  a  magnetic  field  through 
the  use  of  a  fixed  frequency  microwave  source.  At  resonance,  the  magnetoresistance  of  the  Hall  bars  we  made 
showed  an  enhancement  of  the  pxx  component  of  the  resistivity  tensor  (as  expected)  due  to  spin-flips  induced  by 
the  microwave  field.  The  g  factor  was  then  calculated  from  the  field,  B,  and  the  microwave  frequency,  v, 
through  the  simple  formula,  g=hv/pB. 
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Our  second  objective  was  to  demonstrate  the  detection  of  individual  electrons  trapped  electrostatically  in  a 
lithographically  defined  quantum  dot.  Our  approach  was  to  use  a  double-QW  structure  in  which  two 
Ino.53Gao.47As  QWs  surrounded  by  InP  formed  the  confining  layers  for  both  the  electron  qubit  (in  the  upper 
well)  and  the  detecting  layer  (lower  well).  Doping  was  used  above  and  below  the  double-QW  structure  to 
create  a  situation  in  which  carriers  were  localized  primarily  in  the  lower  well.  A  narrow  channel  was  formed  by 
applying  a  negative  bias  to  a  pair  of  depletion  gates,  the  QPC,  coming  in  from  the  side  of  a  mesa  etched  through 
the  double-QW  structure.  Positive  bias  was  applied  to  fonn  the  quantum  dot  via  an  air-bridged  gate  which 
contacted  the  semiconductor  through  a  small  post  typically  less  than  200  mn  in  diameter.  In  the  following 
sections,  we  will  discuss  our  experimental  results  with  respect  to  these  two  objectives. 

(5.2.1)  First  demonstration  of  g  factor  engineered  InGaAs/InP  QWs  with  large  g  factor  shifts 

The  electronic  g  factor  associated  with  electrons 
trapped  in  an  InGaAs  QW  depends  on  the  thickness  of 
the  QW  and  the  composition  of  the  surrounding 
materials  (in  this  case,  InP).  In  the  experiments  we 
have  been  conducting,  30-40  GHz  microwaves  were 
launched  onto  a  microstrip  resonator  fabricated  on  a 
Rogers  5880  laminate  substrate.  The  end  of  the 
microstrip  was  terminated  with  a  short,  and  so,  we 
positioned  the  sample  at  a  location  approximately  one- 
half  wavelength  from  the  end  of  the  microstrip  to 
maximize  the  magnetic  component  of  the  microwave 
field  (Bi)  seen  by  the  sample.  This  configuration 
ensured  the  proper  orientation  of  the  microwave 
magnetic  field  and  the  external  static  field  B0  (namely 
Bj  _L  Bo)  to  observe  electron  spin  resonance  (ESR).  In 
the  experiment,  we  measured  the  component  of  the  g 
tensor  parallel  to  the  static  field  by  monitoring  the 
longitudinal  resistance  of  the  sample  at  resonance.  In 
order  to  observe  the  resonance,  the  sample  must  have  a 
high  mobility  to  allow  spin-splitting  to  be  resolvable  at 
the  resonant  field. 

For  this  experiment,  two  single-QW  samples  were  grown  by  MBE  at  HRL  and  prepared  in  the  Van  der  Pauw 
geometry,  approximately  0.25  in.  square,  in  which  slots  were  cut  using  a  shadow  mask  and  a  sand  blaster.  The 
QW  in  each  of  the  samples  was  2.6  mn  in  thickness  and  composed  of  Ino.44Gao.56As.  The  spacer  thickness  was 

30  mn  for  sample  A  and  40  nm  for  sample  B.  The  resistance  was  measured  using  a  PAR  124 A  lock-in 
amplifier,  with  the  reference  applied  across  the  sample.  For  sample  A,  the  carrier  concentration  was  too  high  to 
observe  spin-splitting  at  an  accessible  field.  Somewhat  serendipitously,  the  carrier  concentration  of  sample  B 
was  just  right,  and  the  resonance  was  observable  in  the  v=\  spin-split  valley  for  microwave  frequencies  between 

3 1  and  40  GHz,  within  the  capabilities  of  our  detection  system. 

For  Fig.  9,  we  measured  the  longitudinal  resistance  (Rxx)  of  sample  B  as  a  function  of  magnetic  field  while  38 
GHz  microwaves  were  introduced  to  the  microstrip  resonator  under  the  sample.  Near  4.4  T,  a  resonant 
enhancement  of  Rxx  was  observed.  The  ESR  signal  (inset  to  Fig.  9)  was  obtained  from  the  raw  data  after 
subtracting  the  non-resonant  background. 


Figure  9:  Longitudinal  resistance  ( Rxx )  vs.  magnetic  field 
for  sample  B  (40  nm  spacer  thickness).  38  GHz 
microwaves  were  applied  using  +10  dBm  source  power, 
producing  a  resonant  enhancement  of  the  resistance 
near  4.4  T.  The  resonant  field  was  obtained  from  the  raw 
data  after  subtracting  the  non-resonant  background 
(inset). 
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In  Fig.  10,  the  ESR  signal  from  sample  B  is  plotted 
for  several  different  microwave  frequencies.  The 
source  power  was  varied  to  maintain  a  sample 
temperature  of  300  mK  during  each  sweep.  ESR 
data  obtained  from  a  phosphorus-doped  Si 
reference  sample  was  used  to  calibrate  the 
magnetic  field.  During  the  calibration  procedure, 
we  measured  a  dependence  of  the  resonant  field  on 
sweep  rate,  due  apparently  to  an  induced  field 
generated  by  eddy  currents  in  the  dewar.  The 
sweep  rate  used  to  obtain  the  data  plotted  in  Fig.  2 
(+0.32  T/min)  was  found  to  generate  an  induced 
field  of  178  G.  The  resonant  field  was  calculated 
by  subtracting  178  G  from  the  peak  position  at 
each  frequency.  The  corrected  values  of  the 
resonance  fields  (Bres)  were  then  used  to  calculate 
the  g  factor  at  each  frequency  according  to  the 
expression  g=h  v/jU/SBres.  The  data  are  in  agreement 
with  the  expression  =  gyo  -  c  (N+Vi)  B, 

where  N  is  the  Landau  level  index.  Our 
measurements  indicate  a  value  of  -0.65  for  g\\0  and 
0.014  1/T  fore. 

The  Overhauser  effect  is  one  possible  source  of 
error  in  our  measurements  of  the  g  factor.  It  is 
particularly  problematic  in  III-V  semiconductors 
since  all  of  the  elements  of  columns  III  and  V  of 
the  periodic  table  have  non-zero  nuclear  spin.  The 
effect  arises  from  an  interaction  between  the 
nuclear  spins  in  the  lattice  and  electron  spins  in  the 
2DEG,  causing  a  temporary  polarization  of  the 
nuclear  lattice  at  resonance.  The  effect  produces 
an  additional  magnetic  field,  Bn,  resulting  in  a  shift 
of  Bres.  The  degree  to  which  the  nuclei  are 
polarized  determines  Bn  and  the  amount  of  shift 
observable  in  Bres.  In  Fig.  1 1,  we  have  plotted  the 
ESR  signal  at  38  GHz  as  a  function  of  microwave 
power.  A  small  shift  (~  200  G)  was  observable  as 
the  power  was  increased  from  1  to  70  mW  at  the 
source.  The  direction  of  the  shift  is  consistent  with 
previous  measurements  of  the  Overhauser  shift  in 
the  AlGaAs/GaAs  materials  system  [4]  although 
its  magnitude  is  much  smaller. 


Figure  10:  Electron-spin  resonance  (ESR)  enhanced 
resistance  plotted  as  a  function  of  field  after  subtracting 
the  non-resonant  background  for  several  applied 
microwave  frequencies.  The  peaks  were  fit  to  a 
Gaussian  function  in  order  to  extract  the  resonant  field. 
An  additional  correction  was  applied  to  account  for  the 
field  induced  during  the  sweep  through  resonance. 
Calculated  g-values  are  shown  in  the  inset  figure. 


Figure  11:  Electron-spin  resonance  (ESR)  enhanced 
resistance  plotted  for  different  source  power  levels  at  38 
GHz.  A  small,  200  G  shift  is  observable,  consistent  in 
direction  with  what  might  be  expected  due  to  the 
Overhauser  effect. 
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(5.2.2)  Fabrication  of  top-gated  air  bridge  device  structures  with  sub-50  nm  feature  sizes 


Our  second  objective  was  to  demonstrate  the  detection 
of  individual  electrons  trapped  electrostatically  under  a 
lithographically  defined  gate.  Our  approach  was  to 
use  the  double  QW  structure  shown  schematically  in 
Fig.  12.  Two  InGaAs  QWs  surrounded  by  InP  form 
the  confining  layers  for  both  the  electron  qubit  (in  the 
upper  well)  and  the  detecting  layer  (lower  well).  Delta 
doping  was  used  above  and  below  the  double  QW 
structure  to  create  carriers  localized  in  the  lower  well, 
resulting  in  a  2DEG.  Input  from  our  modeling  team 
produced  the  nominal  layer  sequence  shown  in  the 
table.  In  the  experiment,  an  applied  bias  was  used  to 
pinch  off  the  detecting  layer  laterally,  nearly  to  the 
point  at  which  it  behaves  like  a  QPC.  The  approach 
was  anticipated  to  make  the  2DEG  highly  sensitive  to 
the  presence  of  individual  charges  nearby.  The 
electron  was  induced  to  tunnel  from  the  detecting 
2DEG  and  into  the  upper  QW  upon  the  application  of 
a  positive  bias  to  an  air-bridged  gate,  deposited 
directly  on  the  AlInAs  cap  layer. 

During  this  course  of  this  program,  we  made 
significant  progress  addressing  some  nagging 
processing  problems  facing  us  in  the  past.  First,  we 
determined  the  annealing  conditions  necessary  for 
creating  low  temperature  ohmic  contacts  to  the  2DEG. 
Second,  we  identified  cleaning  procedures  prior  to 
metal  deposition  that  enhances  adhesion,  making 
possible  contact  to  the  devices  via  wire  bonding. 
Finally,  we  perfected  the  photoresist  layer  deposition 
sequence  and  dosing  parameters  for  electron  beam 
lithography  of  sub- 100  nm  air-bridged  structures  that 
form  the  gate  for  our  single-qubit  logic  devices. 


Alin  As 
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Figure  12:  Schematic  diagram  showing  layer  structure 
and  doping  profile  for  the  nominal  double  quantum  well 
structure  (Q1).  Q2  and  Q3  samples  were  grown 
essentially  to  the  same  parameters  with  2x  and  4x  doping 
levels,  respectively. 


Several  samples  were  grown  by  molecular  beam  epitaxy  according  to  the  layer  sequence  shown  in  Fig.  12.  A 
typical  set  of  three  wafers  were  labeled  Ql,  Q2,  and  Q3.  The  difference  was  in  the  doping  as  Q1  received  the 
nominal  amount,  Q2  received  2x  the  nominal  amount,  and  Q3  received  4x  the  nominal  amount.  One  sample 
(doped  like  Q2)  was  selected  for  dicing  and  subsequent  analysis.  Hall  mobility  and  magnetoresistance  vs. 
magnetic  field  data  confirmed  this  sample  to  be  of  excellent  quality  (mobility  was  >  95,000  cm2/V  s),  exhibiting 
numerous  Shubnikov  de-Haas  oscillations  with  a  carrier  concentration  of  about  4x10  cm'  (about  what  we 
expected).  The  data  are  plotted  in  Fig.  13a,  prior  to  illumination  with  a  red  LED.  Sharp,  spin-resolved  Landau 
level  oscillations  were  observed  at  fields  as  low  as  2.5  Tesla. 


Critical  to  the  performance  of  our  single-qubit  logic  gate  (SQLG)  is  the  requirement  that  only  the  lower  well  be 
filled  below  the  free  surface.  In  order  to  verify  that  this  situation  was  indeed  the  case  for  the  sample  Q2,  we 
illuminated  a  similar  sample  at  low  temperature  and  repeated  the  magnetic  field  sweep.  The  resulting  data  are 
shown  in  Fig.  13b.  Clearly,  the  oscillations  are  considerably  more  complicated.  However,  we  found  that  by 
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Figure  13:  a)  Shubnikov  de-Haas  oscillations  from  a  sample 
doped  like  Q2  exhibit  well-resolved  spin-split  Landau  levels  in  a 
pattern  consistent  with  only  one  well  being  filled  with  electrons 
prior  to  illumination  with  a  red  LED.  b)  After  illumination,  the 
pattern  of  oscillations  changes  as  additional  carriers  become 
available  to  fill  the  well  closer  to  the  sample  surface.  A 
straightforward  Fourier  analysis  of  the  data  plotted  vs.  inverse 
field  is  inset  and  shows  a  second  periodicity  corresponding  to 
the  electron  density  in  the  upper  well  (peak  at  1.8X1011  cm"2). 


Figure  14:  a)  Model  of  the  HRL  single-qubit  logic 
gate  (SQLG)  showing  quantum  point  contact  gates 
(L  and  R)  and  the  gate  electrode.  The  qubit  is 
represented  by  an  electron  drawn  into  the  upper 
quantum  well  by  the  quantum  gate.  The  presence 
of  the  electron  is  detected  through  the  influence  of 
its  charge  on  the  conductance  of  the  lower 
quantum  well  in  the  region  constricted  by  the 
quantum  point  contact  gates,  b)  SEM  micrograph 
of  a  typical  device  from  the  wafer  Q1  shows  space 
under  the  air  bridge  and  slight  undercutting  of  the 
quantum  point  contact  gates.  The  gate  diameter  is 
approximately  200  nm. 


taking  a  Fourier  transform  of  the  data  plotted  as  a  function  of  inverse  field,  two  periodicities  are  observed.  The 
second  periodicity  arises  from  the  fact  that  upon  illumination,  additional  carriers  are  created  and  trapped  in  the 
upper  well.  Prior  to  illumination,  only  one  periodicity  is  observed,  hence,  we  were  confident  that  the  samples 
were  grown  to  the  proper  specifications  and  the  doping  spanned  the  range  of  values  we  intended.  Ql,  Q2,  and 
Q3  were  then  submitted  to  the  process  engineers  for  fabrication  into  SQLG  devices. 


Prototype  SQLG  devices  were  fabricated  on  long  rectangular  mesas  etched  down  to  the  InP  substrate  onto 
which  a  thin  (30  nm)  SfiNx  dielectric  was  deposited.  Ohmic  contacts  were  annealed  at  either  end  of  the  mesa, 
along  with  side  contacts  used  to  eliminate  lead  resistances  from  the  measurements.  A  layer  of  Ti/Pt/Au  was 
deposited  onto  and  into  a  hole  etched  in  the  dielectric  at  the  center  of  the  mesa.  Tri-level  photoresist  was 
deposited  on  the  entire  structure  and  electron-beam  lithography  was  used  to  develop  a  mold  into  which  more 
Ti/Pt/Au  metallization  was  deposited  to  form  the  air-bridged  gate  structure.  A  close-up  view  of  the  gate  area  is 
shown  in  Fig.  14.  In  Fig.  14a,  a  model  showing  the  configuration  of  the  QPC  gates  (L  and  R)  and  the  quantum 
dot  gate  electrode  (air  bridge  structure)  can  be  seen.  Current  flows  right  to  left  in  the  figure  through  the  QPC 
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gates  under  normal  operation.  Negative  bias  is  applied  to  the  QPC  gates,  while  positive  bias  on  the  quantum 
gate  draws  the  electron  qubit  into  the  upper  well.  Fig.  14b  shows  a  scanning  electron  micrograph  (SEM)  image, 
taken  from  one  of  the  devices  fabricated  from  Ql.  The  dimension  of  the  gate  is  roughly  200  nm.  The  smallest 
gates  we  made  measured  about  50  nm  in  diameter.  A  slight  undercutting  of  the  QPC  gates  is  also  seen  in  Fig. 
14b.  A  modification  to  the  electron-beam  lithography  exposure  sequence  was  found  during  the  processing  of 
Q2  to  virtually  eliminate  this  problem. 


(5.2.3)  Demonstration  of  single-electron  sensitivity  i 

As  long  as  the  electrons  occupying  the  quantum  dot 
remain  in  place  long  enough  (on  the  order  of  several 
hundred  milliseconds)  our  detection  setup  should  be 
capable  of  detecting  single-electron  tunneling  events. 
For  this  reason,  we  increased  the  spacer  thickness  in 
the  nominal  structures  from  the  10.6  nm  used 
previously  in  Ql,  Q2,  and  Q3  to  30  nm  in  the  present 
design.  The  expectation  was  that  as  electrons 
tunneled  into  the  dot  sequentially,  we  would  observe 
discontinuous  jumps  in  the  conductance  of  the 
channel  due  to  Coulomb  repulsion.  This  is  exactly 
what  we  have  now  observed. 

In  Fig.  15,  we  plotted  differential  conductance 
(dl/dVsD)  as  a  function  of  gate  voltage  for  a  50  nm 
quantum  dot  located  at  the  center  of  a  1  pm  wide 
QPC.  The  data  were  taken  using  a  PAR  124 A  lock- 
in  amplifier  with  the  reference  applied  through  a 
resistor  divider  network.  In  this  way,  22  nA  AC 
current  was  applied  to  the  sample  at  190  Hz.  With  - 
3.85  V  applied  to  the  QPC  and  zero  bias  on  the  gate 
(dot),  this  particular  device  was  pinched  off  (zero 
conductance),  as  can  be  seen  in  the  figure.  By 
increasing  the  gate  voltage,  the  channel  became 
conducting.  Above  0.8  V,  a  noticeable  change  in  the  si 


top-gated  InGaAs/InP  quantum  dot  devices 


Figure  15:  Differential  conductance  (dl/dVSo)  as  a 
function  of  gate  voltage  for  a  50  nm  dot  in  a  1  pm  gap 
between  two  leads  of  a  quantum  point  contact  (QPC). 
The  data  were  obtained  for  a  double-quantum  well 
structure  containing  InGaAs  wells  in  an  InP  matrix. 
Reproducible  structure  (plateaus)  is  notable  in  addition  to 
“noise”  appearing  at  0.8  V  and  higher  in  gate  voltage. 
The  QPC  was  held  at  a  large  negative  bias  to  deplete  the 
channel  under  the  dot. 

;nal-to-noise  ratio  was  observed. 


By  slowing  down  the  sweep  rate  and  taking  many  more  points,  the  “noise”  observed  in  the  data  for  Fig.  15 
began  to  look  much  more  like  steps.  In  Fig.  16,  we  show  four  sweeps  over  the  same  range  in  gate  voltage  with 
the  QPC  voltage  still  fixed  at  -3.85  V.  The  data  reveal  discontinuous  jumps  to  lower  values  of  differential 
conductance.  These  jumps  can  be  identified  with  single-electron  transitions.  Since  few  events  were  observed 
below  0.7  V  gate  bias,  we  believe  that  the  quantum  dot  was  indeed  in  the  few-electron  occupancy  regime. 


In  Fig.  16(b),  we  have  expanded  the  scale  to  show  more  clearly  the  individual  levels  corresponding  to  the  zero 
and  higher  electron  states  of  the  quantum  dot.  The  transitions  are  clearly  separated  in  gate  voltage  and 
reproducible,  from  trace-to-trace.  Impressively,  one  can  follow  the  transitions  on  any  particular  trace  and  count 
electrons  tunneling  into  the  dot  one-by-one,  up  to  as  many  as  15  electrons!  This  is  quite  remarkable  and 
demonstrates  a  level  of  control  that  we  were  extremely  excited  to  see. 


16 


Figure  16:  (a)  dl/dVSD  obtained  from  the  same  device 
as  described  in  Fig.  15  at  a  significantly  reduced  sweep 
rate  and  increased  point  density.  The  data  reveal 
transitions  between  multiple  levels  of  conductance 
which  we  tentatively  associate  in  (b)  with  the  zero 
electron,  one  electron,  and  two  electron  states  of  the 
quantum  dot.  With  the  arrival  of  each  successive 
electron,  the  conductance  of  the  quantum  well  below  it 
decreases  in  discrete  steps. 


Figure  17:  Differential  conductance  vs.  time  under  fixed 
bias  conditions  shows  random  fluctuations  in  the  dot 
occupancy  by  ±  one  electron  (random  telegraph  signal). 
The  average  time  between  events  is  on  the  order  of 
several  minutes,  due  to  the  rather  large  (30  nm)  spacing 
between  the  dot  and  the  lower  quantum  well  in  this 
particular  device  structure.  Two  multiple-electron  events 
are  also  observed  and  are  now  believed  to  be  caused  by 
unintentional  spikes  in  the  gate  voltage. 


In  another  set  of  experiments,  we  measured  (in  a 
somewhat  qualitative  way)  the  tunneling  rate  for 
electrons  as  the  gate  voltage  was  swept  through  the 
various  transitions.  In  one  test,  we  allowed  the  device  to 
reach  equilibrium  at  some  (high)  value  of  gate  voltage 
for  which  several  electrons  were  allowed  to  tunnel  into 
the  dot.  We  then  quickly  reduced  the  gate  voltage  to 
<0.7  V,  where  the  dot  should  have  emptied.  Instead,  we 
observed  the  conductance  to  rise  in  steps  over  a  period 
of  several  minutes,  until  finally,  all  the  electrons  had 
tunneled  out  of  the  dot.  To  confirm  that  the  electrons 


were  gone,  we  further  lowered  the  gate  voltage  and 
waited  at  least  20  minutes  without  seeing  any  more 
events.  These  tests  have  convinced  us  that  we  were  dealing  with  extremely  long  tunneling  times  (on  the  order 
of  several  minutes)  and  there  were  only  a  few  electrons  on  the  dot. 


In  Fig.  17,  we  plotted  differential  conductance  as  a  function  of  time  for  fixed  gate  voltage  and  QPC  voltage. 
The  gate  was  held  at  +1.0  V,  a  rather  high  value  for  which  several  (perhaps  7  or  8  electrons)  were  on  the  dot.  In 
the  figure,  one  can  see  the  conductance  switch  back  and  forth  between  basically  two  levels  in  a  random  pattern 
reminiscent  of  a  random  telegraph  signal  (RTS)  but  over  a  long  timescale.  This  behavior  can  be  attributed  to 
single  electrons  tunneling  into  and  out  of  the  dot  once  every  few  minutes  or  so.  The  two  spikes  appear  to  be 
events  during  which  suddenly  all  of  the  electrons  are  forced  to  jump  out  of  the  dot  and  return  one-by-one  to  the 
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previous  level  of  occupation.  These  events  have  been  traced  to  momentary  electromagnetic  interference  or  line 
voltage  fluctuations  which  caused  100  ps  long  spikes  in  the  gate  voltage. 

(5.2.4)  First  demonstration  of  few-electron  quantum  dot  spectroscopy  using  random  telegraph  analysis 

During  the  last  year  of  the  program,  we  focused  our 
attention  on  the  single  event  transitions  we  had 
observed  previously  in  our  SQLG  devices.  Recall  that 
in  these  devices,  electrons  are  induced  to  tunnel 
through  an  InP  spacer  layer  from  the  detecting  2DEG 
and  into  the  quantum  dot.  Therefore,  decreasing  the 
spacer  thickness  should  increase  the  tunneling  rates. 

For  the  devices  made  using  the  30  nm  InP  spacer,  the 
tunneling  rates  were  too  long,  making  a  quantitative 
study  nearly  impossible.  Here,  we  describe 
experiments  on  a  device  made  from  a  structure  with 
only  a  15  nm  thick  spacer,  resulting  in  much  faster 
tunneling  rates  (one  event  every  second  or  less). 

The  faster  tunneling  rates  in  these  devices  means  that 
some  transitions  may  not  appear  as  discrete  jumps  but 
rather  as  small  inflections  in  the  channel  current  as  a 
function  of  gate  voltage.  Discontinuous  jumps  are 
only  resolvable  if  the  measurement  bandwidth  is 
higher  than  the  lifetimes  of  the  particular  states  of  the 
quantum  dot  involved  in  the  transition.  In  our 
experiments,  we  used  the  reference  signal  from  a  Stanford  Research  SR830  digital  lock-in  amplifier  to  source  a 
small  AC  current  at  2  kHz  through  the  channel  under  our  quantum  dot.  We  used  a  voltage  divider  to  provide  a 
40  pVrms  AC  excitation  to  one  end  of  the  channel  and  connected  the  other  end  to  the  current  input  on  the  lock- 
in  (virtual  ground).  The  lock-in  time  constant  was  set  to  30  ms  and  therefore  limited  the  bandwidth  of  the 
measurement  to  a  little  more  than  30  Hz.  Signal-to-noise  ratio  considerations  prevented  us  from  increasing  the 
bandwidth  any  further. 

In  Fig.  18,  we  show  the  current  as  a  function  of  gate  voltage  for  a  SQFG  device  with  a  1  pm  wide  QPC  and  a 
100  nm  diameter  dot  electrode.  The  voltage  on  the  QPC  was  set  to  -4.30  V.  Clearly  visible  are  inflections  and 
occasional  discrete  jumps  associated  with  single-electron  tunneling  from  the  channel  into  the  quantum  dot. 
The  slowest  transitions  are  shown  in  the  three  insets  with  the  transition  near  1.2  V  corresponding  to  an  N  to 
N+l  transition  with  N~8.  This  particular  transition  was  slow  enough  for  us  to  resolve  individual  tunneling 
events  as  we  stepped  the  gate  voltage  and  sampled  the  current  as  a  function  of  time  to  generate  a  random 
telegraph  signal  (RTS).  Due  to  the  thinner  spacer  layer,  timescales  are  now  short  enough  to  gather  good 
statistics  and  extract  quantitative  information.  The  signal  to  noise  ratio  was  sufficient  to  allow  for  extraction  of 
lifetime  data. 

In  Fig.  19,  we  plotted  the  current  as  a  function  of  time  for  three  gate  voltages  spanning  the  width  of  the 
transition  near  1.2  V.  In  each  case,  the  current  was  observed  to  jump  between  two  levels,  corresponding  to  a 
transition  from  the  N  electron  state  (higher  current)  to  the  N+l  electron  state  (lower  current)  and  vice-versa. 
The  data  for  Vgate  =  1.174  V  are  offset  slightly  for  clarity.  The  full  set  of  data  were  obtained  by  stepping  the 
gate  voltage  in  2  mV  steps  between  1.160  and  1.180  V  and  sampling  the  current  in  1  sec  intervals  for  1  hour. 
At  Vgate  =  1.170  V,  the  lifetime  of  the  N  and  N+l  electron  states  are  roughly  equal.  For  voltages  less  than 


Gate  Voltage  (V) 


Figure  18:  Current  vs.  gate  voltage  for  a  SQLG  device 
fabricated  using  a  structure  with  a  15  nm  spacer  between 
the  detecting  2DEG  and  the  quantum  dot.  The  QPC 
voltage  was  set  to  -4.30V.  The  data  were  taken  by 
applying  the  reference  signal  of  an  SR830  lock-in 
amplifier  through  a  voltage  divider  to  provide  an  AC 
excitation  of  40  pVrms.  The  frequency  and  the  lock-in 
time  constant  were  set  to  2  kHz  and  30  ms,  respectively. 
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1.170  V,  the  dot  is  primarily  in  the  N  electron  state 
and  for  voltages  greater  than  1.170  V,  it  is  mostly 
in  the  N+l  electron  state.  The  distribution  of 
lifetimes  was  determined  to  obey  simple  Poisson 
statistics.  This  is  a  well  understood  phenomenon, 
observed  for  years  in  MOSFET  traps  [5]  and  more 
recently  in  quantum  dots  [6]. 

An  analysis  of  the  data  set  allowed  us  to  extract  a 
variety  of  information  about  the  transitions  and 
further  supports  our  conclusion  that  electrons  are 
tunneling  between  the  gate-induced  quantum  dot 
and  the  readout  channel  below.  We  began  by 
writing  an  algorithm  to  ignore  the  relatively  small 
background  drift  in  the  current  and  output  a  list  of 
the  length  of  time  the  dot  is  in  one  state  vs.  the 
other.  The  distribution  of  times  were  binned  and 
plotted  in  terms  of  the  probability  that  the  state  is 
N  or  N+l  for  a  given  length  of  time.  On  a  semilog 
plot,  the  data  for  each  gate  voltage  are  roughly 
linear,  as  expected  for  Poisson  statistics  and  the  fit 
to  the  data  provides  a  measure  of  the  lifetime  of  that 
ratio. 


Random  Telegraph  Signals  for  V05-010, 15  nm  spacer 
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Figure  19:  Current  as  a  function  of  time  at  fixed  gate 
voltage.  At  Vqate  =  1.170  V,  the  lifetimes  for  the  N  electron 
state  and  the  N+1  electron  state  of  the  quantum  dot  are 
roughly  equal.  For  VGate  =  1.164  V  (1.174  V),  the  number 
of  electrons  on  the  dot  is  primarily  N  (N+1).  The  distribution 
of  lifetimes  obeys  simple  Poisson  statistics. 

state.  The  lifetimes  are  plotted  in  Fig.  20  along  with  their 


As  a  function  of  gate  voltage,  we  found  that  only  the  ratio  of  the  lifetimes  varies  monotonically.  The  lifetimes 
actually  increase  anomalously,  near  the  middle  of  the  transition.  We  surmised  that  this  phenomenon  reflects 
kinetics  as  opposed  to  thermodynamics.  In  equilibrium,  the  ratio  for  each  gate  voltage  should  be  predictable 
and  reflect  the  relative  occupancy  of  the  states.  In  fact,  the  occupancy  probability  is  easily  derived  from  a  grand 
canonical  ensemble  formalism  Z  =  where  P=l/kT,  EF  is  the  Fermi  energy  (chemical  potential), 

{s,N} 

gs  is  the  degeneracy  (orbital  and/or  spin)  of  the  state  s,  and  the  summation  is  over  all  possible  electronic  states 

and  possible  numbers  of  electrons,  N.  Then,  the  probability  of  configuration  I  is  p(Ni,Ei )  =  -^gle  IHErN‘EF) . 

Now  consider  a  system  with  two  states  of  charge  involving  energy  levels  N  and  N+l.  Assume  En  =  0 
(arbitrary)  and  take  En+i  to  be  the  energy  at  which  the  system  transitions  to  the  N+l  state.  Then  the  partition 
function  for  this  two  state  system  is  Z  =  gNet>NF‘F  +  gN+xe  ,I(Fn"  (A"' Then  the  probability  of  occupying  the 


state  N+l  is  given  (with  g  =  )  by  pN+i  =  [1  +  gep(-Emx  £,r1"1  1 


8 


]  .  The  ratio  of  lifetimes  can  then  be  interpreted 


JV+l 


1  Pn+ i 


as  the  occupancy  ratio  —  =  - — ±'N+l  =  g  exp(( EQD  -  EF)/ kT .  This  ratio  is  the  exponential  in  the  energy 

xi  Pn+ i 

difference  between  the  N+l  quantum  dot  state  and  the  Fermi  energy.  The  fit  to  the  experimental  data  assumes 
that  the  energy  difference  is  linear  in  the  gate  voltage  over  this  narrow  range. 


Increasing  the  gate  voltage  moves  the  energy  level  of  the  dot  relative  to  the  Fermi  energy.  A  fit  to  the  data  in 
Fig.  20  tells  us  how  the  energy  level  of  the  dot  state  relative  to  the  Fermi  level  varies  with  gate  voltage.  This 
can  be  compared  with  device  simulation  results.  We  measure  therefore  a  AEqd(Vgate)  =  9  meV  *  Vgate  at  300 
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mK.  This  compares  to  the  theoretical  prediction  of 
12  meV  *  Vgate-  This  is  completely  equivalent  to 
measuring  the  width  of  a  Coulomb  blockade  peak. 


We  also  have  calculated  tunneling  rates  associated 
with  the  single  electron  charging  events  in  these 
devices.  Tunneling  rates  were  computed  using 
Fermi's  Golden  Rule  based  on  numerical  eigenstates 
derived  from  fully  3 -dimensional  self-consistent 
Poisson-Schroedinger  calculations.  The  effects  of 
Coulomb  interaction  on  the  quantum  dot  states  and 
hence  tunneling  rates  was  considered.  A  comparison 
of  these  numerical  results  with  experimental 
estimates  of  tunneling  rates  derived  from  random 
telegraph  signals  indicates  that  the  observed 
tunneling  rates  are  considerably  slower  than 
expected  from  theory.  Detennining  the  cause  of  this 
discrepancy  will  be  part  of  a  future  research 
program. 

In  summary,  RTS  has  been  observed  in  enhancement 
mode  quantum  dots  formed  by  a  single  gate  placed  in  forward  bias.  All  data  are  consistent  with  single  electron 
charging  of  a  quantum  dot.  This  demonstrates  basic  proof  of  principle  of  this  device.  The  RTS  rate  for  some  of 
the  transitions  is  in  the  range  that  allows  for  good  statistics  to  be  obtained.  The  tunneling  rates  are  slow  enough 
to  be  resolved  by  our  several  tens  of  Hz  bandwidth  measurement  system,  and  fast  enough  to  allow  for 
accumulation  of  good  statistics  in  a  reasonable  amount  of  time  (hours).  Analysis  of  the  data  shows  that  the 
lifetimes  are  Poisson  distributed,  as  expected.  The  ratio  of  occupancies  was  found  to  be  exponential  in  gate 
voltage  and  the  width  of  the  transition  (inferred  from  the  slope)  agrees  with  theoretical  estimates  assuming  a 
300  mK  electron  temperature. 


Lifetimes  of  High  and  Low  Current  States 


Figure  20:  Lifetime  of  the  N  electron  state  (high  current) 
vs.  the  N+1  electron  state  (low  current)  as  a  function  of 
gate  voltage  for  the  transition  studied  in  Fig.  19.  A 
monotonic  decrease  in  the  ratio  of  the  lifetimes  is  observed, 
with  a  slope  consistent  with  device  simulation  results. 


(5.3)  Control  accomplishments 

The  initial  control  effort  was  centered  on  developing  a  simulation  model  of  the  essential  features  of  the 
proposed  spin-coherent  photon  transmitter/receiver  system.  The  purpose  of  the  model  is  to  evaluate  control 
strategies  for  creating  the  requisite  quantum  logic  gates.  In  particular,  the  control  design  goal  is  to  make  an 
EPR  pair  -  the  Bell  states.  Input  preparation  and  photon/electron  interactions  were  not  initially  included. 
Emphasis  was  placed  on  modeling  sources  of  uncertainties  in  the  system,  e.g.,  uncertain  parameters,  fringe 
effects,  and  sources  of  decoherence  such  as  nuclear  spin. 

The  control  variables  in  the  model  are  the  external  potentials  (gate  voltages)  which  directly  effect  the  g-factors 
in  the  semiconductor  material  in  the  presence  of  an  external  magnetic  field.  The  g-factors  are  used  to  control  the 
electron  spin  states.  Preliminary  results  indicated  that  some  implementations  of  the  logic  gates  needed  to 
execute  the  teleportation  algorithm  place  constraints  on  the  design  of  the  device. 

Various  control  strategies  were  evaluated  in  simulation  in  order  to  map  out  a  robust  design  space.  The  design 
approaches  included  pulse  methods  similar  to  those  used  in  NMR  as  well  as  optimal  control  design.  The  latter 
approach  is  much  easier  to  implement  in  this  device  than  in  NMR,  for  example,  because  the  driving  fields  are 
applied  to  gate  voltages  rather  than  magnetic  fields. 
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In  the  second  year  of  the  program,  we  examined  the  electron/nuclei  interaction  as  a  source  of  decoherence. 
Proposed  models  found  in  the  literature  produce  decoherence  through  an  interaction  that  behaves  like  a 
measurement,  i.e.,  the  electron/nuclei  interaction  is  dependent  on  the  magnitude  of  the  electron  wave  function. 
Assuming  this  interaction  is  relatively  small,  we  studied  the  use  of  linearized  models  to  obtain  optimal 
corrections  to  control  pulses  designed  without  the  interaction. 

We  also  concentrated  on  developing  and  refining  methods  for  quantum  tomography  and  and  quantum  state 
detection  A  number  of  problems  in  quantum  state  estimation  (state  tomography),  quantum  system  identification 
(process  tomography)  and  quantum  state  and  system  detection  can  be  cast  as  convex  optimization  problems. 
The  great  advantage  of  convex  optimization  is  a  globally  optimal  solution  can  be  found  efficiently  and  reliably, 
and  perhaps  most  importantly,  can  be  computed  to  within  any  desired  accuracy  using  an  interior-point  method. 

Some  of  these  problems  are  already  known  to  be  convex  but  have  not  fully  exploited  the  available  convex 
solvers  or  duality  theory.  For  example,  it  is  known  that  Maximum  likelihood  Estimation  (MLE)  of  the  quantum 
state  (density)  is  a  convex  optimization.  What  we  also  showed  is  that  a  number  of  other  MLE  problems  are 
convex,  e.g.,  estimating  the  distribution  of  known  states  and  quantum  process  tomography  in  the  Kraus  operator 
sum  representation  (OSR)  in  a  fixed  basis.  One  important  problem  that  is  not  convex  is  MLE  of  Hamiltonian 
parameters.  We  showed,  however,  how  duality  theory  can  help  establish  bounds  on  the  parameter  estimates  for 
this  problem. 

Another  problem  that  can  be  solved  via  convex  optimization  is  experiment  design.  (Experiment  design  here 
means  choosing  the  number  of  experiments  to  be  perfonned  in  a  particular  system  configuration;  a 
configuration  being  any  number  of  combinations  of  sample  times,  hardware  settings,  etc.  For  example,  in 
quantum  state  photonic  tomography,  we  can  determine  the  optimum  number  of  wave  plate  settings  to  achieve  a 
desired  estimation  accuracy.)  We  have  been  able  to  apply  the  experiment  design  procedure  invoked  by  the 
Cramer-Rao  Inequality  to  all  the  MLE  problems  mentioned  above,  including  MLE  of  Hamiltonian  parameters. 
In  all  these  cases,  the  optimum  experiment  design  problem,  although  integer-combinatoric,  can  be  relaxed  to  a 
convex  optimization  problem  whose  solution  provides  upper  and  lower  bounds  on  the  unknown  optimal  integer 
solution.  The  MLE  of  the  state  or  process  can  be  combined  with  the  optimal  experiment  design  in  a 
"bootstrapping"  iteration  to  make  the  estimation  more  efficient. 

We  have  addressed  the  problem  of  designing  a  detector  that  is  maximally  sensitive  to  specific  quantum  states. 
The  design  problem  can  be  formulated  as  a  convex  optimization  problem  in  the  matrices  of  the  POVM  (positive 
operator  valued  measure)  that  characterize  the  measurement  apparatus,  or  with  a  given  POVM,  the  matrices 
which  characterize  the  OSR  in  a  fixed  basis.  We  specifically  addressed  maximizing  the  posterior  probability  of 
detection  and  showed  that  this  is  a  quasiconvex  optimization  problem  in  either  the  POVM  or  OSR  matrices. 
Previous  work  in  this  area  has  only  considered  the  joint  probability  of  detection  over  POVM  matrices, 

In  all  the  cases  described  above  we  showed  how  duality  theory  can  be  used  in  various  special  cases  to  give 
insight  into  the  nature  (and  possible  physical  implementation)  of  the  optimal  solutions. 

To  make  quantum  information  systems  robust  will  require  that  the  systems  operate  under  some  form  or 
combination  of  control  and  error  correction.  Moreover,  it  is  highly  likely  that  in  order  to  achieve  the  desired 
system  objectives,  these  systems  will  have  to  be  tuned,  or  even  entirely  detennined,  using  estimated  quantities 
obtained  from  data  from  the  actual  system  rather  than  solely  relying  on  an  initial  design  from  a  theoretical 
model.  To  alleviate  these  concerns,  we  have  been  developing  the  mathematical  and  computational  tools  for 
quantum  state  and  process  tomography  as  well  as  quantum  state  detection. 
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For  quantum  state  detection  the  goal  is  to  design  the  measurement  apparatus  to  determine  whether  or  not  a 
particular  state  was  present  in  the  input  to  the  detector.  The  problem  is  thus  equivalent  to  the  design  of  a 
quantum  channel  that  optimally  transforms  the  input  state  distribution  (assumed  to  be  given,  and  including  non- 
orthogonal  states)  to  the  output  measurement  outcomes.  The  channel  may  be  lossy,  and  may  introduce  noise, 
and  thus  there  may  be  latency  in  the  measurement,  in  which  certain  outcomes  are  ambiguous.  Previous  work  in 
this  problem  has  shown  that  maximizing  the  average  conditional  probability  over  the  set  of  all  possible 
measurement  devices  is  a  convex  optimization,  specifically  what  is  called  an  SDP,  or  Semidefinite  Program. 
We  have  found  that  maximizing  the  posteriori  probability  over  all  possible  measurement  apparatuses  is  also 
convex,  specifically  quasiconvex  and  therefore  also  solved  via  convex  programming.  The  posteriori 
probability,  which  answers  the  question:  if  detector  declares  that  state  is  present,  what  is  the  probability  that  it  is 
actually  present.  We  developed  software  for  solving  this  problem. 

We  also  investigated  the  problem  of  generating  the  control  signals  to  activate  the  quantum  logic  gates.  At 
present  we  see  this  as  a  classical  control  and  signal  processing  problem.  The  main  difficulty  is  that  the 
interconnects  become  very  dense  as  they  approach  the  geographical  quantum  gate  location  and  this  density 
gives  rise  to  “cross-talk”  across  the  different  channels.  More  precisely,  at  the  gigahertz  frequencies  required  for 
quantum  gate  operation,  the  interconnect  system  behaves  exactly  like  a  multi-input  multi-output  transmission 
line.  The  problem  is  further  compounded  by  uncertainties  arising  from  manufacturing  variations  in  the  quantum 
gate  impedances  as  well  as  in  the  interconnects. 

We  are  specifically  interested  in  determining  the  limit  of  performance  possible  for  such  a  system.  We  explored 
several  approaches: 

(1)  Open-loop  robust  control 

Design  the  controls  using  prior  knowledge  or  estimates  of  the  interconnect/quantum  gate  system  and  uncertainty 
ranges. 

(2)  Open-loop  adaptive  control 

Design  the  controls  using  knowledge  or  estimates  of  the  interconnect/quantum  gate  system  obtained  from  on¬ 
line  measurements  and  experiments  specifically  created  for  obtaining  such  knowledge. 

(3)  Closed-loop  control 

Use  real-time  feedback  from  the  interconnect/quantum  gate  system  to  regulate  the  control  signals  to  the  gates. 

To  make  the  problem  more  concrete  we  starting  with  the  assumption  that  the  interconnect  system,  not  including 
the  quantum  gate  impedances,  can  be  very  well  modeled  or  estimated  from  on-line  data.  The  main  source  of 
uncertainty  is  thus  the  impedance  to  the  quantum  logic  gates.  Even  knowing  the  interconnect  system,  however, 
does  not  mean  signals  can  be  created  to  perfectly  invert  the  dynamics.  There  is  a  limit  to  bandwidth  and  signal 
timing  that  is  possible  and  certainly  made  more  difficult  when  the  quantum  gate  impedences  are  uncertain.  One 
saving  grace  is  that  the  gate  uncertainties  are  likely  to  be  systematic  rather  than  continuously  random  variations. 
Thus  once  identified,  the  problem  becomes  one  of  inversion. 

Towards  this  end  we  concentrated  on  approach  using  closed-loop  control.  The  first  issue  is  how  to  generate  a 
feedback  signal.  Adding  more  interconnects  to  obtain  measurements  obviously  compounds  the  problem  in 
several  ways.  First  it  increases  the  possibility  of  more  uncertainty  and  more  cross-talk.  Secondly,  it  requires 
more  geography,  which  may  reduce  cross-talk.  A  simpler  approach  is  to  take  advantage  of  the  transmission  line 
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bi-directional  capabilities.  In  this  regard  we  have  conceived  of  a  method  to  use  the  reflected  waves  as  feedback. 
This  has  the  disadvantage  of  adding  signal  delay.  However,  the  system  uncertainty  has  not  been  increased 
thereby  making  robust  inversion  more  possible,  modulo  the  additional  delay.  We  need  also  account  for 
uncertainties  from  this  measunnent  scheme  arising  from  noise  and/or  uncertain  measunnent  dynamics.  The  goal 
of  this  investigation  was  to  not  only  understand  the  performance  limits  but  also  to  develop  design  tools  which 
can  be  utilized  for  any  proposed  system. 

Finally,  we  have  applied  methods  of  convex  optimization  to  the  design  of  quantum  error  correction  (QEC). 
Convex  optimization  has  proven  very  powerful  in  our  previous  efforts  in  state  and  process  tomography  [quant- 
ph/041 1093]  and  quantum  state  detection  [quant-ph/0403 150].  For  QEC  it  remains  an  open  question  as  to 
whether  errors  can  be  dealt  with  more  efficiently  in  tailored  system  configurations  that  do  not  require  the  full 
level  of  complexity  attendant  to  standard  QEC. 

We  have  shown  that  the  problem  of  designing  a  quantum  infonnation  error  correcting  procedure  can  be  cast  as  a 
bi-convex  optimization  problem,  iterating  between  encoding  and  recovery,  each  being  a  semidefinite  program 
(SDP)  [quant-ph/0606078].  For  a  given  encoding  operator,  the  problem  is  convex  in  the  recovery  operator.  For 
a  given  method  of  recovery,  the  problem  is  convex  in  the  encoding  scheme.  This  allows  us  to  derive  new  codes 
that  are  locally  optimal.  We  now  have  examples  of  such  codes  that  can  handle  errors  that  are  too  strong  for 
codes  derived  by  analogy  to  classical  error  correction  techniques.  We  have  also  shown  that  the  dual 
optimization,  also  an  SDP,  is  of  lower  complexity  and  thus  requires  less  computational  effort.  The  SDP 
formalism  also  allows  for  a  robust  design  by  enumerating  constraints  associated  with  different  error  models. 

A  distance  measure  was  also  developed  [quant-ph/0606064]  between  two  unitary  propagators  of  quantum 
systems  of  differing  dimensions  along  with  a  corresponding  method  of  computation.  A  typical  application  is  to 
compare  the  propagator  of  the  actual  (real)  process  with  the  propagator  of  the  desired  (ideal)  process;  the  former 
being  of  a  higher  dimension  then  the  latter.  This  measure  has  the  advantage  of  dealing  with  possibly  correlated 
inputs,  but  at  the  expense  of  working  on  the  whole  space  and  not  just  the  information  bearing  part  as  is  usually 
the  case,  i.e.,  no  partial  trace  operation  is  explicitly  involved.  This  distance  measure  and  an  average  measure  of 
channel  fidelity  both  depend  on  the  size  of  the  same  matrix:  as  the  matrix  size  increases,  distance  decreases  and 
fidelity  increases. 
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